Anomalous behavior of spin wave resonances in Gai_ x Mn x As thin films 
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We report ferromagnetic and spin wave resonance absorption measurements on high quality epi- 
taxially grown Gai-^Mn^ As thin films. We find that these films exhibit robust ferromagnetic 
long-range order, based on the fact that up to seven resonances are detected at low temperatures, 
and the resonance structure survives to temperatures close to the ferromagnetic transition. On 
the other hand, we observe a spin wave dispersion which is linear in mode number, in qualitative 
contrast with the quadratic dispersion expected for homogeneous samples. We perform a detailed 
numerical analysis of the experimental data and provide analytical calculations to demonstrate that 
such a linear dispersion is incompatible with uniform magnetic parameters. Our theoretical analysis 
of the ferromagnetic resonance data, combined with the knowledge that strain-induced anisotropy 
is definitely present in these films, suggests that a spatially dependent magnetic anisotropy is the 
most likely reason behind the anomalous behavior observed. 

PACS numbers: 75.50.Pp, 76.50.+g, 75.30.Ds, 75.30.Gw 



I. INTRODUCTION 

There has been great deal of interest in elucidating the 
properties of Gai-^Mn^As and other diluted magnetic 
semiconductors (DMSs) during the past several years. 
The intense research on these materials is partly mo- 
tivated by the fact that they hold promise as building 
blocks of "spintronic" semiconductor devices^ 2 * 3 -. In- 
deed, the incorporation of magnetic properties in semi- 
conductor heterostructures can, in principle, lead to the 
development of new devices that manipulate both the 
spin and the charge degrees of freedom of the carriers. 
Due to its compatibility with conventional electronic de- 
vices, Gai-zMn^As is one of the most readily usable alloy 
systems for exploring spintronic prototypes. However, 
these materials exhibit a series of fascinating strong cor- 
relation phenomena that are not yet fully understood, 
such as a metal-insulator transition, field-induced ferro- 
magnetism, and magnetoresistance 3 -. Before any appli- 
cations become possible, it is necessary to provide a de- 
tailed description of both electronic and magnetic struc- 
ture and properties of this class of materials. Investigat- 
ing the magnetic properties of the ferromagnetic ground 
state in Gai-^Mn^As, with Curie temperatures as high 
as T c ~ 160^, is therefore especially important. 

Given the fact that DMSs arc synthesized in film form 
by using molecular-beam epitaxy (MBE), ferromagnetic 
resonance spectroscopy (FMR) is the most suitable ex- 
perimental probe for studying the dynamics of the ferro- 
magnetic order parameter, that also allows for the spec- 
troscopy of the spin wave excitations. FMR is a powerful 
technique to study magnetic properties in magnetic thin 
films^ and has already been used by several groups in the 
magnetic characterization of the Gai-^Mn^As films££. 
The same technique can be used to obtain the resonance 
fields of the spin wave modes in a thin film. Extract- 
ing the various magnetic parameters influencing the spin 
excitations is essential for gaining complete control of 



the magnetic properties of DMS films with an eye on 
successful future applications, such as spin injection and 
manipulation. 

The main objective of this paper is to provide a self- 
consistent picture of the properties of spin wave excita- 
tions based on a comparison of experimental results and 
theoretical calculations. The theory of spin wave reso- 
nance (SWR) has been developed four decades ago by 
Kittefi, who also pointed out that FMR measurements 
are capable of detecting several modes of the magnetic 
excitations in ferromagnetic thin films. For an external 
magnetic field normal to the surface of a homogeneous 
film he finds the following resonance condition for spin 
waves pinned at its boundaries. 
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where H n is the external field value (at fixed external 
radiation frequency) where the n-th SWR mode is ob- 
served, Ho is the magnetic field that corresponds to the 
ferromagnetic resonance, L is the sample thickness and 
the constant D is proportional to the stiffness constant 
defined in Sec. EE- 

In this paper we report the observation of such spin 
wave resonances in Gai-^Mn^As. These experiments 
provide a direct proof for true long-ranged ferromag- 
netic order in Gai-^-Mn^As. Surprisingly, the spin waves 
which we observe exhibit a somewhat unusual behavior: 
We find a spin wave spectrum with H n ~ n. A re- 
cent work by Goennenwein and collaborators^ reported 
H n ~ n 2 / 3 . Both results arc in qualitative disagreement 
with Eq. QJ, which states that for a homogeneous film 
the resonance field H n of the n th mode should be propor- 



tional to 



We trace back the origin of the anomalous 



dispersion to the magnetic properties of Gai-^Mn^As 
thin films, and find that in order to understand the 
anomalous spin wave dispersion, it is crucial to allow the 
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magnetic parameters inside the film to depend on the 
distance from the interface. In principle, such inhomo- 
geneity in the profile of magnetization, spin stiffness, or 
magnetic anisotropy^ could all result in such an anomaly. 
However, our experimental and theoretical results pre- 
sented in this paper point towards the presence of uni- 
axial anisotropy and/or spin stiffness that depends on 
the distance z from the surface of the film. More specifi- 
cally, we find that our resonance experiments can be well 
understood by assuming a uniaxial magnetic anisotropy 
with a quadratic dependence on the distance z from the 
Gai-^Mn^As/GaAs interface . The presence of strong 
magnetic anisotropy in Gai-^Mn^As and Ini-^Mn^As 
epitaxial films has already been clearly demonstrated by 
a series of experiments^*S*ifi'ii*iSiii, which also indicate 
that the magnetic anisotropy depends on the lattice mis- 
match between the substrate and the DMS layer. 

This paper is organized as follows: In Section ^ we 
present the experimental results for a group of as-grown 
and annealed Gai„ x Mn ;r As samples and discuss briefly 
the position, intensity and linewidth of the observed res- 
onance peaks. In Section II I II we present a general spin 
wave equation of motion that allows for the variation 
of the magnetization and magnetic anisotropy along the 
film thickness, and discuss our numerical calculation for 
the position and intensity of the resonance peaks. Fi- 
nally, in Section ll VI we compare the experimental results 
with the theoretical calculations and discuss the possi- 
ble explanations for the inhomogeneity of the magnetic 
anisotropy. 



II. EXPERIMENTAL RESULTS 

Recently, a systematic study of the fundamental FMR 
mode has been reported for a series of Gai-^Mn^As films 
grown on various substrates^. The dependence of the 
FMR position on the angle between the applied magnetic 
field and the crystallographic axes of the sample was care- 
fully documented, and detailed information has been ob- 
tained on the magnetic anisotropy and its variation with 
temperature. The uniaxial and cubic anisotropies deter- 
mined experimentally generally corroborate with earlier 
theoretical predictionsii*i£. In this paper the same exper- 
imental technique is used to study spin wave resonances 
(SWRs) in thin Gai_ x Mn x As films. 

Although DMSs have a very low concentration of mag- 
netic atoms, and these system are often described in 
terms of percolation^ and impurity band models^*^, 
it was possible to observe SWR spectra with up to seven 
spin wave mode a 9 ' 17 , demonstrating that real long-range 
order develops in these materials. 

The three Gai-^Mn^As thin films analyzed in this pa- 
per have previously been studied by Sasaki et al, who 
observed both the SWR and the FMR lines. Their ini- 
tial work concentrated on the uniform mode (the FMR 
line), and they used this feature to investigate the over- 
all magnetic anisotropy in Gai-^Mn^As films4. All the 
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TABLE I: Experimental parameters extracted from, and 
used in the theoretically analysis of, the spectra of three 
Gao.924Mno.076 As samples of different thicknesses, pre- and 
post-annealing. H max is the resonance field of the highest 
(ferromagnetic resonance) peak, M is the static magnetiza- 
tion measured by SQUID, T c is the critical temperature, and 
AH is the linewidth. When the width varies with the mode 
number (see text), the range of linewidths is given. 

samples were grown by low temperature molecular beam 
epitaxy (LT-MBE) on semi- insulating GaAs substrates. 
The Mn concentration x = 0.076 was determined by x- 
ray diffraction, and the Curie temperature and depen- 
dence of magnetization dependence on applied magnetic 
field and temperature were obtained by superconduct- 
ing quantum interference device (SQUID) measurements. 
The FMR and SWR measurements were carried out us- 
ing a 9.46 GHz microwave spectrometer. 

The dc magnetic field H ext was oriented perpendicular 
to the film plane (see Fig. QJ. As a result of the lattice 
mismatch between the substrate and the Gai-^Mn^As 
films, the samples exhibit compressive strain in the sam- 
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FIG. 1: Schematic diagram of the experimental geometry. 
The Gai_iMn x As film is grown on a thick GaAs substrate. 
A large constant magnetic field H ex t is applied in the z direc- 
tion: H ex t\\z. This field produces a large magnetization M(z) 
which points predominantly in the z direction. Additionally, 
a small external microwave magnetic field h _L H ext produces 
a small harmonically-varying perturbation in magnetization, 
m _L z 
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d= 200 nm 
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FIG. 2: (a) FMR and SWR absorption spectra for as-grown 
Gao.924Mno.o76As samples of different thickness (100 nm, 
150 nm and 200 nm) at T=4 K. (b) FMR and SWR spectra 
for the 200 nm thick Gao.924Mno.o76As sample for different 
temperatures, before and after annealing. 



pie plane, which leads to a strong uniaxial anisotropy. 
Consequently, in the absence of the external field the 
magnetization lies in the plane parallel to the film. How- 
ever, the applied static field at which FMR and SWRs 
are observed is strong enough to align the magnetization 
perpendicular to the sample surface. After measuring 
the SWR in the as-grown samples, we annealed them 
for 60 minutes. As a result, the Curie temperature T c 
was increased by about ~ 40% and the magnetization 
by ~ 25%. Table [|] summarizes some of the charac- 
teristic properties of three representative Gai-^Mn^As 
samples. The values listed in the table were extracted 
from SQUID, FMR and SWR measurements. Clearly, 
the annealing process has a profound effect not only on 
the critical temperature and magnetization of the films, 
but on practically all characteristics of the the ferromag- 
netic and spin wave resonances, including the tempera- 
ture and thickness dependence of the resonance position 
and linewidth. In order to obtain more insight into the 
behavior of these materials, we now proceed to evaluate 
the actual FMR spectra before and after annealing. 

Fig. [3 shows typical examples of the observed spectra. 
The main trends and qualitative features of these data 
can be effectively conveyed if the position and linewidth 
of the resonance features are plotted as a function of the 




FIG. 3: Experimentally observed T = 4K SWR mode po- 
sitions (a), linewidths (b) and peak intensities (c) for the 
200 nm thick Gao.924Mno.076 As sample before (open circles) 
and after (open squares) annealing as a function of mode in- 
dex n. (a) the solid lines show the linear fitting for the ex- 
perimental data, (c) The dotted line connects the theoretical 
result for the normalized peak intensities. Note that this cal- 
culation was done with no fitting parameters (see text). Odd 
modes are not observed experimentally, as their amplitudes 
are much smaller compared to those of the even modes. 



mode number. We label the largest resonance at the high- 
est field (the uniform FMR mode) as n = 0, and the SWR 
modes as n = 1,2, 3..., where the increasing mode index 
corresponds to decreasing resonance field. The results of 
this procedure are shown in Fig. [3] The distance between 
SWR modes increases for thinner samples. This is in ac- 
cordance with the simple picture, that these resonances 
arc standing waves of magnetization trapped between 
the two interfaces of a uniform thin film, and that the 
level spacing between successive standing wave modes in- 
creases with decreasing geometrical size of the 'resonator 
cavity'. However, from this simplistic picture also fol- 
lows that the resonances must depend quadratically on 
the mode index, in accordance with Kittcl's equation, QJ. 
This expectation fails spectacularly in our Gai-^Mn^As 
thin films. As one can see in Fig. [3] for the 200 nm thick 
as-grown and annealed samples, in spite of the fact that 
we are able to see as many as 7 spin- wave resonances, the 
SWR mode positions as a function of the mode number 
do not follow the expected quadratic law. Instead, the 
resonance fields H n exhibit a linear dependence on the 
mode index n. 

We can also compare the intensity of the resonances 
with Kittel's predictions. However, to do that we have 
to keep in mind that the FMR measurement is done with 
a field modulation lock-in technique, which measures the 
derivative of the absorption as a function of the applied 
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magnetic field. The observed peak-to peak height AI n 
and peak-to-peak width AH n of the n-th resonance are 
related to its intensity I n aa*& 

I n ~ AI„ AH 2 n (2) 

It is rather striking that the lincwidth behaves qual- 
itatively differently for annealed and unannealed sam- 
ples: For the as-grown film the linewidth is indepen- 
dent of n and its value is AH n ~ 200 Oe. For the an- 
nealed samples, on the other hand, AH n decreases with 
n (see Fig. 01b), which is a typical behavior of metallic 
thin filmsi^. This change to more conventional behav- 
ior in the resonance lincwidth corroborates the observa- 
tion that the annealed samples are more metallic than 
the as-grown samples: resistivity decreases upon anneal- 
ing, consistent with an increase in the density of mobile 
charge carriersSS. On the other hand, the peak intensities 
remain anomalous despite the annealing. The extracted 
intensities /„ are shown in Fig.|3Jc. They follow a similar 
law for both the annealed and unannealed samples, and 
the decrease of the successive intensities is much slower 
than the /„ oc l/(n+l) 2 dependence predicted by Kittcl. 

Clearly, based on the qualitative and quantitative ex- 
perimental observations made above, the Gai-^Mn^As 
films do not behave as homogeneous ferromagnetic thin 
films. The first attempt to explain this behavior is 
due to Sasaki et alJUL who, inspired by the work of 
Portis22, suggested a picture where the magnetization 
inside the ferromagnetic thin film is not homogeneous in 
the direction perpendicular to the plane of the film. Al- 
though such inhomogeneity could explain qualitatively 
the linear dependence of successive SWR positions on 
n, the gradient of composition required to fit the data 
is unrcalistically large, and recent neutron reflectivity 
measurements 2 ^ do not support this degree of variation of 
magnetization across the sample. An alternative model 
is therefore needed for comparing experiment and the- 
ory. We present such a theoretical analysis in the fol- 
lowing section, which will retain in spirit the approach 
of Refs. Il7ll2ll in that we also propose inhomogeneity in 



In Eq. ©, M = |M| denotes the magnetization of the 
sample at temperature T, and H ext \\z is the applied dc 
magnetic field. Every integration runs over the volume V 
of the film. Note that the free energy functional in Eq. © 
only describes transverse fluctuations of the magnetiza- 
tion, and the longitudinal fluctuations are assumed to be 
negligible, M = |M| = constant. 

The first term in Eq. @ is an exchange free energy, 
with A(r) denoting the spin wave stiffness, while the sec- 



magnetic parameters along the thickness of the sample 
as the main reason behind the observed anomalous be- 
havior of the SWR positions. However, in our theoretical 
framework - in addition to a gradient of magnetization 
- we also allow for a variation in magnetic anisotropy 
and spin stiffness. One of the attractive features of our 
approach is that - once the anisotropy profile has been 
determined to fit the resonance positions - it allows for 
a parameter-free evaluation of the normalized peak in- 
tensities (see the points connected by dashed lines in 
Fig. E|c. )• The results are in excellent agreement with 
those extracted from the experiment for both as-grown 
and annealed samples 2 ^. This agreement in particular 
gives us confidence that the theoretical model and numer- 
ical analysis which we present below gives an essentially 
accurate description of the FMR and SWR experiments 
in Gai_ x Mn x As. 

III. THEORETICAL APPROACH 

A. Semiclassical spin wave equations 

It is well accepted that the magnetism of the 
Gai-zMujAs is due to indirect interaction between Mn 
spins mediated by holesi^2Si2Si2i. Although this system 
contains considerable positional disorder - which, com- 
bined with spin-orbit effects, may lead to non-collincar 
ground states2L2£ - it appears that on larger length 
scales, relevant for the long wavelength collective modes 
of the ferromagnetic order parameter, these DMS materi- 
als behave as conventional ferromagnets^. In the follow- 
ing we will therefore neglect many of the complications 
which are not relevant when dealing with spin waves, 
and we shall employ the usual semi-classical equations of 
motion to study spin wave excitations in the absence of 
damping2tt. 

The first step to derive the semiclassical equations of 
motion is to construct the free energy functional for the 
magnetization M(r). For temperatures below the tran- 
sition temperature Tq this can be expressed as 



(3) 

I 

ond term represents - depending on the sign of A'(r) - 
a uniaxial or an in-plane anisotropy energy with respect 
to direction u. The primary source of the anisotropy K 
is the strain field due to the lattice mismatch between 
the GaAs substrate and the ferromagnetic Gai-^Mn^As 
filmSi. Since the Gai_ x Mn x As films discussed through- 
out this paper were deposited on GaAs, they have an 
in-plane easy axis which corresponds to u = (0, 0, 1) and 
positive K. We have not included other cubic anisotropy 



J 



F(M) = d 3 r A | VM| 2 + J d 3 r A |M • u | 2 - J d 3 r (H ext - 2ttM) ■ M . 
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terms in Eq. because at temperatures they are small 
compared to \K\. Finally, the last term in Eq. [21 is the 
Zeeman free energy together with the demagnetization 
energy. Note that all coefficients in Eq. Q depend on 
temperature T, and on spatial coordinates. 

Given the material specifics detailed above, the semi- 
classical equations of motion can be written as follows: 



dM(f, t) 
dt 



| 7 |M(r,t) xH t ot(r,t). 



(4) 



Here 7 is the gyromagnctic ratio, which we approximate 
by that of the Mn core spins, 7 « -f e — 2/xb/?i. H tot (r, t) 
denotes the total effective magnetic field, which can be 
obtained by taking the functional derivative of the free 
energy functional, Eq. with respect to the magneti- 
zation, H tot (r,<) = -6F({M})/5M(r,t). The result at 
T = K temperature is 



H tot (r,t) = -AV 2 M(r,t) 



2K 
~M~' 



2A 
TP 
first 



4ttM, (5) k x ,k y 



where we assumed that the exchange constant A = 
does not depend on the spatial coordinate. The 
term in Eq. 03 is usually called the exchange field, while 
the second term defines the uniaxial anisotropy field, 



a 2K '« 



To compute the spin wave spectrum for a field parallel 
to z we expand the magnetization around its equilib- 
rium value and orientation Mz and, as we mentioned be- 
fore, we allow for transverse deviations only: M = Mz + 
m(r, t). Assuming periodic time dependence and plane 
wave character of the r.f. magnetization in the x and y 
directions, m(r, t) = [?7i x (z)x + m y (z)y] e tuJt e' lk:cX+ ' lkyy , 
one obtains the following equation of motion for 
m + (z) = m x (z) + irriy(z) in the long wavelength limit 



0: 



{—AM 



dz 2 



AwM(z)- H a {z) 



. d 2 M{z), . 
A 2 ; }m + (z) = (H ext 

I 



-)m + (z) , 



(6) 



where w r is the frequency of the microwave radiation. 
In course of the derivation of Eq. HJ1 we assumed that 
only the magnetic anisotropy and the magnetization vary 
along the Z direction (growth direction) (H a = H a (z) 
and M = M(z)), but (as already mentioned after Eq. 
(5)) that A(z) = A is constant. While these assump- 
tions need not be valid in general, we will show below 
that the experimental data on our DMS films can be 
well explained by the presence of a spatially-dependent 
anisotropy fields. As mentioned before, a spatially vary- 
ing magnetization can also reproduce qualitatively the 
main features of the datai^, but if one allows only M(z) 
to vary, we now find by numerically solving Eq. an 
unrealistic solution that cannot reproduce quantitatively 
the experimental SWR spectrum. There is yet another 
possibility, that the exchange constant is itself spatially 
dependent, A = A(z). Given that A depends on the 
carrier concentration, a scenario with a z-dependent ex- 
change constant would, in principle, be consistent with 
recent experiments of Koeder and coworkers^. While 
we cannot exclude this possibility, we expect that this 
would give a spectrum similar to that of the spatially 
non-uniform magnetization, and therefore a very large 
spatial dependence of the carrier concentration would be 
needed to explain the experimental data. We therefore 
classify this possibility, together with that of spatially- 
dependent magnetization, as sub-dominant mechanisms, 
which would at most play a secondary role in explaining 
the features of our resonance experiments. The versatil- 



ity of our numerical analysis is clearly evident in evaluat- 
ing the different scenarios discussed above: our numerical 
scheme allows us to check all these scenarios against our 
resonance absorption data and the constrains set by other 
experiments, and select the most viable one for explain- 
ing our FMR and SWR measurements. 

Equation JBJl has the form of a one-dimensional 
Schrodinger equation for an electron in a quantum 
well. If we compare the coefficients in EqEl with the 
Schrodinger equation, we can see that H 2 /(2AM(z)) is 
the analog of the electron mass (in our case it may de- 
pend on the position z if M varies with position) and 



4irM(z)-H a (z)-A 



d'M(z) 
dz'i 



is the analog of the potential 



energy V(z). Equation © can then be formally written 

as 



[^+V(z)]mt(z)=E n m+(z), (7) 

where p = —iti-^ is the momentum operator and E„ = 
H ex t — — is the n-th eigenvalue. Having solved this 
Schrodinger equation, we can compute the intensities I n 
of each mode as2i 



\Jo m t( z )dz\ 2 
Jo \mt{z)\ 2 dz' 



(8) 



where the integration runs over the ferromagnetic film 
thickness. 
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B. Boundary conditions and intensity of SWR 

To solve Eq. ©, it is crucial to establish the bound- 
ary conditions for SWR. While different boundary con- 
ditions do not affect the positions of the resonances in an 
essential way, the intensity of the SWR peaks strongly 
depends on them. It is therefore important to choose the 
appropriate boundary conditions (BCs). 

The question of boundary conditions has been debated 
for a rather long time in the literature and, to our knowl- 
edge, no general agreement has been reached. Different 
authors22ii&Si treat the BC problem in different ways, 
and assume different boundary conditions appropriate to 
the specific sample which they are studying. 

Equation iJBJ must be, in general, solved together with 
macroscopic Maxwell equations. Ament and Rado^i ar- 
gued that in the absence of any magnetic anisotropy the 
spin wave solutions should satisfy so-called free (or anti- 
node) boundary conditions { dm d }^ = at z = and 
z = L). On the another hand, Pincus^ pointed out 
that Ament and Rado have not treated properly the dis- 
continuity of the space at the surface of the film. Pin- 
cus (and also previously Kitted) indicated that at the 
interface an additional term, proportional to the gradi- 
ent of the magnetization at the surface (not included in 
equation |BJ), is also allowed by symmetry, and that this 
term can lead to the pinning of spin waves at the bound- 
aries. Furthermore, Pincus observed that free BCs can 
never be appropriate when there is surface anisotropy. 
He also showed that an antiferromagnetic oxide layer on 
the surface of a film can give rise to a surface anisotropy 
which pins the spins at the end-points (m£(z = 0) = 
and 171+ (z = L) = 0). Finally, Pincus and Kittel have 
shown that it is more appropriate to use the so-called 
dynamic boundary conditions'^ which for the lowest 
lying modes effectively reduce to pinned boundary con- 
ditions. 

We believe that in the case of Gai_ x Mn x As films one 
must also use pinned boundary conditions. This is be- 
cause, first, there is a strong strain field present in these 
films that generates a large anisotropy field, and sec- 
ond, we expect a strong surface-induced anisotropy in 
the vicinity of the surface due to the strong spin-orbit 
coupling in Gai_ x Mn x As, which would also pin the spin 
waves«i. Fortunately, our numerical analysis allows us 
to try different boundary conditions. We were unable 
to obtain a good fit to the peak intensities when we 
used unpinned or partially unpinned boundary condi- 
tions. Pinned boundary conditions, on the other hand, 
gave very good agreement with the measured SWR spec- 
tra. This is illustrated in Fig. 01 where we show our 
best fits with pinned and unpinned boundary conditions 
for the SWR spectra for an as-grown 200 nm sample at 
T = 4 K. 

Typical solutions for pinned boundary conditions and 
the corresponding intensities are shown in Fig. [S] The 
intensity of odd modes is suppressed, because of the can- 
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FIG. 4: The derivative of the absorption of the SWR as a 
function of applied dc magnetic field for an as-grown 200 nm 
Gai_ I Mn I As film. The solid lines represent the experimental 
data, and the dashed lines show our theoretical results for the 
magnetic anisotropy profile of Eq. @. The upper/lower pan- 
els show the results of our computations for pinned/unpinned 
boundary conditions. 

cellation of the regions with m + > and m + < 0, re- 
spectively. In fact, in Fig. 0] the odd modes are barely 
visible, and only peaks associated with the even modes 
can be observed. 

The low energy spin (standing) waves tend to be local- 
ized around the region with higher magnetic anisotropy. 
Qualitatively speaking, this is because we have the nega- 
tive of the spatially-dependent uniaxial anisotropy act- 
ing as the effective trapping potential for the modes, 
and therefore it is easier to create spin waves where the 
anisotropy is large. For smaller values of A, spin waves 
tend to be more localized around the large anisotropy 
points. On the other hand, spin waves with higher en- 
ergies will be extended over the whole sample thickness, 
and their energy will not depend linearly on n. 



IV. NUMERICAL RESULTS AND DISCUSSION 

We solved Eq. (jSJ) numerically, and varied the exchange 
parameter A and the profiles of the anisotropy field H a (z) 
(or the magnetization M(z)) to obtain a best fit to the 
experimental spectra. From these calculations we con- 
cluded that to explain a linear variation of the resonance 
fields, H n ~ n, the changes of H a (z) (or M(z)) must be 
substantial across the whole depth of the sample. For 
profiles where the variation of H a (z) (or M(z)) is only 
near the surface, the resonance positions show always a 
quadratic behavior, H n ~ n 2 . This can be easily under- 
stood from Eq. J7J), since in this case the potential V(z) 
is similar to that of a rectangular quantum well and the 
eigen-modes are therefore quadratically separated. 

As discussed in the previous section, the spatial varia- 
tion of the magnetization and the exchange constant are 
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FIG. 5: Ilustration of the functional form of three spin wave 
modes with n = 0, 1 and 2 (solid, dashed, and dotted lines) 
for the as-grown L — 200 nm sample at T — 4K. We used an 
exchange constant A = 3000 nm 2 and a = 0.8. The rescaled 
anisotropy profile is marked by open circles. Distances are 
measured in units of L. In the inset we show the calculated 
intensities I n of the modes in the main figure . 

sub-dominant processes. On the other hand, we were able 
to obtain quantitative agreement with the experimental 
data allowing the magnetic anisotropy field H a = H a (z) 
to vary across the sample while keeping the magnetiza- 
tion constant, M{z) = const. Therefore, in the discus- 
sion that follows, we shall mostly focus to the case of a 
z-dependent anisotropy field, H a = iJ a (z), and assume 
M = const. 

The computed SWR spectrum depends on the specific 
shape of the anisotropy profile H a (z). We found, in par- 
ticular, that a linear dependence of H a on z, H a (z) ~ z, 
is clearly in disagreement with our experimental data. 
However, we could obtain an excellent fit by assuming a 
quadratic dependence on z, H a {z) ~ z 2 . More specifi- 
cally, both of the following profiles fit the experimental 
data rather well: 

H a (z) = ff Q cxp (-az 2 /L 2 ) (9) 
= H a (l-az 2 /L 2 ), , (10) 

where L is the thickness of the film, a ~ 0.75 is a di- 
mensionless fitting parameter, and the maximum value 
of the anisotropy field, H a (0), has been extracted from 
the analysis of the main FMR resonance^. Since we can 
obtain the value of the saturation magnetization M from 
accurate SQUID measurements, we end up with only two 
fitting parameters, a and the stiffness A (or equivalcntly, 
the exchange constant A). 

The origin of the linear behavior on n can be easily un- 
derstood from the formal analogy with the Schrodingcr 
equation (Eq. EJ) and the energy spectrum of the har- 
monic oscillator: as long as the difference between 



the resonance field H n and the maximum value of the 
anisotropy field |.ff a (0)| — H n is small compared to 
we expect the corresponding wave function m + (z) to be 
well approximated by the Hcrmite functions and to have 
a linear behavior ~ n. This condition is well satisfied for 
the experimentally observed spin wave resonance fields. 

In order to study the anisotropy profile, we focused 
on the 200 nm sample, because this sample exhibited 
the largest number of resonance peaks. First, we chose 
our parameters to fit the 4 K data (Fig. |HJ): Having 
obtained the magnetization M ~ 17.9 emu/cm 3 from 
SQUID experiments, and the magnetic anisotropy field 
H a = 4.400 Oe from the angular dependence of the main 
FAIR resonance line, we adjusted A = 3000 ± 300 mn 2 
and a = 0.75 to obtain the best fit which corresponds 
to a spin stiffness A = 0.4 meV/A. For the spec- 
trum measured at 40 K in the as-grown sample we ob- 
tained a very good fit with the same exchange con- 
stant A and a, but replacing H a (0) by the measured 
40 K anisotropy H^ 0K (0) = 4.200 Oe and magnetization 
M i0K = 9.4 emu/cm 3 . In both cases the SWR spectra 
(and thus the resonance fields H n ) exhibit a quasi-linear 
dependence on n. 

Annealing of Gai-^Mn^As has various - presently 
not entirely understood - effects. First, the satura- 
tion magnetization increases upon annealing. Further- 
more, Rutherford backscattering experiments suggest^ 
that the primary process resulting from annealing is the 
diffusion of interstitial Mn ions out of the film. Since 
interstitial Mn takes away carriers from the hole band, 
annealing results probably in an increase of the carrier 
density, and the resulting carrier density is most likely 
inhomogeneous. Since both the anisotropy and the ex- 
change energy are sensitive to the carrier density, an- 
nealing is expected to change the value of both H a and 
A^i and it may also change their profiles. 

Consequently, it is natural to expect that the SWR 
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FIG. 6: Comparison between theoretical (solid lines) and 
experimental (symbols) data for the 200 nm sample. Two 
different temperatures for the as-grown sample and one tem- 
perature for the annealed sample are considered. The satura- 
tion magnetization data from SQUID measurements is used 
in the theoretical approach. 
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spectra observed on annealed specimens will require fit- 
ting parameters that are different from the ones used for 
the as-grown sample. Nevertheless, we were able to ob- 
tain a good fit to the experimental SWR by keeping the 
profile unchanged (a = 0.75), adjusting the anisotropy 
field i/* nnoal (0) = 5300 Oe and using A = 1225 nm 2 . Sur- 
prisingly, this value of exchange constant A corresponds 
to the same value of the stiffness parameter A as in the 
as-grown sample. 



V. CONCLUSIONS 

We have analyzed the experimental FMR and SWR 
spectra of thin Gai-^Mn^As films grown on GaAs sub- 
strates by low temperature MBE technique. We com- 
pared the experimentally observed resonance positions 
and intensities with theoretical calculations. The exper- 
imental results are clearly inconsistent with the assump- 
tion of a homogeneous Gai-^Mn^As film. After solving 
the semiclassical equations of motion and using indepen- 
dent experimental data for the magnetic anisotropy and 
magnetization, we have shown that a nearly quadratic z- 
dependence of the uniaxial anisotropy, H a (z) — H a (fi) ~ 
z 2 could adequately explain the position and intensity of 
the experimental data. An inhomogencous magnetiza- 
tion could also lead to similar effects. However, the spa- 
tial variation and the magnitude of magnetization needed 
to explain the experimental data is incompatible with the 
experiments. Spatial variation of the magnetization, as 
observed in Ref. play a secondary role in the SWR 
spectra. From our numerical analysis we also concluded 
that the observed spin waves are pinned at the surfaces 
of the film. As a consequence of pinning, the intensity of 
n = odd spin wave modes is suppressed and the experi- 
mentally observed spectra consist of even modes only. 

Goennenwein et al. performed similar SWR mea- 
surements and also found spin wave spectra that are 
incompatible with having homogeneous Gai-^MnajAs 
filmsS. They, however, find that H n oc n 2 / 3 , consis- 
tent with a linear z-dcpcndcncc of the anisotropy field, 
H a {z) — H a (0) ~ z. This is in qualitative contrast 
with our findings for the resonance fields which scale as 
H n ~ n. Furthermore, Goennenwein et al. appear to 
have implicitly assumed free boundary conditions, which 
implies that both even and odd modes are observable. In 
contrast, we argue that pinned boundary conditions are 
more adequate for our samples, and we find that, again, 
only even modes can be observed. 

There may be several mechanisms that produce 
inhomogencous magnetic parameters. The uniaxial 
anisotropy is primarily due to uniaxial strain field in the 
film that develops due to the lattice mismatch between 
the substrate and the ferromagnetic filmS*i2iii. In prin- 
ciple, an inhomogeneous strain field could therefore pro- 
duce an inhomogeneous anisotropy field H a . However, in 
MBE grown samples the strain usually relaxes suddenly, 
when the sample reaches a critical thickness, where dislo- 



cations start to form and relax the strain. In these non- 
equilibrium MBE grown samples, however, there seems 
to be no strain field relaxation at all through dislocation 
formation: Even for micron thick samples the measured 
in-plane lattice constant of the film is the same as that 
of the substrate. 

One can also obtain inhomogencous magnetic proper- 
ties by assuming an inhomogeneous hole concentration. 
Indeed, annealing experiments support the notion that 
the dominant compensation effect is due to interstitial 
Mn ions, which act as double donors and may also com- 
pensate the spin of substitutional Mn. These intersti- 
tial Mn ions can diffuse out of the sample during the 
growth process, which takes usually a few hours and 
takes place at the same temperature as the annealing. 
As a result, it is quite possible that the concentration 
of charge carriers (related to that of Mn intcrstitials) 
varies across the film. Since the exchange anisotropy 
and exchange constants are both related to the carrier 
density^, this can serve as a mechanism to produce a 
z-dependent anisotropy field. Several recent experiments 
provide firm support for an inhomogeneous hole concen- 
tration: Koeder and coworkers^ find indications of gra- 
dients in both the carrier concentration and Curie tem- 
perature of epitaxial Gai-^Mn^As films. The presence 
of such gradients are also in accord with recent obser- 
vations of interstitial Mn diffusionSS^. The notion that 
Mn diffusion is a key process for annealing-induced en- 
hancement of magnetism has recently been backed up by 
annealing experiments on samples capped^ with a few 
monolayers of undoped GaAs, which would block the in- 
terstitial Mn from diffusing out of the sample. Under 
these circumstances the annealing process gave no sig- 
nificant increase in the Curie temperature. We believe 
that these experimental results, together with our reso- 
nance measurements and the theoretical arguments given 
above provide solid support for our calculation that a 
gradient in the magnetic parameters must be present in 
Gai^Mn^As. 

We also studied the linewidth of the observed reso- 
nances. For the as-grown samples the SWR linewidths 
do not depend on the mode number. This behavior hints 
that the relaxation is due to spin-orbit coupling^ In- 
deed, the measured (rather large) resonance width is 
compatible with the presence of a relatively large ran- 
dom anisotropy. 

For the annealed samples, on the other hand, the 
linewidth decreases with increasing mode number, which 
is characteristic to eddy current relaxation in metallic 
samples 19 . This behavior is consistent with the results 
of annealing experiments, since for higher carrier densi- 
ties the random anisotropy effects become less important, 
and at the same time the sample becomes more metal- 
lic. It also underscores the fact that there is a significant 
physical difference between as-grown and annealed sam- 
ples, going beyond quantitative changes in the values of 
saturation magnetization and the Curie temperature. 

Since the variation of clastic and/or magnetic proper- 
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ties across the Cai^Mn^ As film can have important con- 
sequences in its future spintronics applications, a more 
detailed experimental and theoretical analysis is neces- 
sary to understand and control the magnetic properties of 
these thin ferromagnetic layers. One way to obtain more 
information about the effects of the surface anisotropy is 
to perform FMR measurements in symmetric three-layer 
structures GaAs/Gai-aMn^As/GaAs and compare them 
with the previously obtained results. A systematic study 
of the thickness and annealing time dependence of sim- 
ilarly grown samples would be also important to under- 
stand the origin of the observed gradient of composition. 
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